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The Extended Theory of Finite Fermi Systems is based on the conventional Landau- Migdal the- 
ory and includes the coupling to the low-lying phonons in a consistent way. The phonons give rise 
to a fragmentation of the single-particle strength and to a compression of the single-particle spec- 
trum. Both effects are crucial for a quantitative understanding of nuclear structure properties. We 
demonstrate the effects on the electric dipole states in 208 Pb (which possesses 50% more neutrons 
then protons) where we calculated the low-lying non-collective spectrum as well as the high-lying 
collective resonances. Below 8 MeV, where one expects the so called isovector pygmy resonances, 
we also find a strong admixture of isoscalar strength that comes from the coupling to the high-lying 
isoscalar electric dipole resonance, which we obtain at about 22 MeV. The transition density of this 
resonance is very similar to the breathing mode, which we also calculated. We shall show that the 
extended theory is the correct approach for self-consistent calculations, where one starts with effec- 
tive Lagrangians and effective Hamiltonians, respectively, if one wishes to describe simultaneously 
collective and non-collective properties of the nuclear spectrum. In all cases for which experimental 
data exist the agreement with the present theory results is good. 
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I. INTRODUCTION 

The structure of neutron-rich nuclei is important for 
many applications in astrophysics. The rapid neutron- 
capture process, for instance, requires detailed knowledge 
of the properties of nuclei between the valley of stabil- 
ity and the neutron drip-line, such as neutron capture 
rates or photon-induced neutron emission cross sections. 
Progress in Nuclear Resonance Fluorescence photon scat- 
tering experiments has made possible studies of energet- 
ically low lying strength [l|, 0, JljEl and to distinguish 
between 1~, 2+ and 1 + states Surprisingly, in the 

same region, appreciable isoscalar electric dipole strength 
was already detected some years ago with (a, a' 70) ex- 
periments and reported in [6[. New data using the same 
technique have already been analyzed 

In analogy to the giant dipole resonance which ex- 
hausts the major part of the Thomas- Reiche-Kuhn sum 
rule, the low lying strength is called the pygmy resonance. 
Commonly, this strength is interpreted as an oscillation 
of the neutron halo against the nuclear core Q . A collec- 
tive model predicts the pygmy B(E1) strength to increase 
with the neutron excess [9(. 

For the electromagnetic dissociation of light nuclei 
with a one-neutron halo, Typel and Baur have devel- 
oped a model independent approach that relies on only 
a few low-energy parameters, e.g. the neutron sep- 
aration energy. Coulomb dissociation is found to be 



dominated by the energetically low-lying electromagnetic 
dipole strength. While in experimentally known nuclei, 
the relevant low-energy parameters can be obtained di- 
rectly from the data, an extrapolation to unknown nuclei 
remains a challenge for many-body theory [Iol |. 

Several microscopic models have been applied to the 
pygmy resonances. The self-consistent calculations are 
based on an effective Lagrangian or Hamiltonian which 
has the advantage that the extrapolation from the known 
stable nuclei to the nuclei near the drip line is well- 
defined. Theories of this kind are quite successful in 
reproducing the nuclear masses and radii, and also the 
collective excitation modes of nuclei, such as the giant 
dipole resonance, using the quasi-particle random phase 
approximation. A recent review on the relativistic ap- 
proach is given in Ref.[ll[. Quantitatively, the relativis- 
tic approaches to the pygmy resonances do not compare 
too well with the experimental data. The theory ob- 
tains mean excitation energies of the pygmy strength 
consistently above the experimental strength. This ap- 
pears to be a general phenomenon observed by several 
groups. Goriely and Khan [l^] calculated within the 
QRPA the El strength distributions for all nuclei with 
8 < Z < 1 10 between the proton and neutron drip lines 
using known Skyrme forces. In their calculation the low- 
lying El strength was located systematically higher by 
some 3 MeV compared with the available data. 

Here, we want to point out that the discrepancy be- 
tween the theoretical and experimental pygmy strength 
may be overcome by pushing the many-body approach 
beyond the mean field approximation. 

In the past twenty years theoretical investigations of 
the structure of nuclei are performed in self-consistent 
models as mentioned before and in Migdal's theory of fi- 
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nite Fermi system (TFFS). In both cases one ends up 
with the random phase approximation(RPA), which de- 
scribes the excited states of the nuclei as superpositions 
of particle-hole states. The input for the RPA equations 
are single-particle energies, single-particle wave functions 
and the residual particle-hole (ph) interaction. The two 
approaches differ in the way how the input data are deter- 
mined. In the self-consistent approach one starts with an 
effective Lagrangian or effective Hamiltonian and calcu- 
lates the single-particle energies and single-particle wave 
functions in the mean field approximation. Also, the 
residual interaction is derived from the original interac- 
tion. This approach is considered the most fundamental 
one, as the effective Lagrangians and effective Hamilto- 
nians can be used for all nuclei. As most of the self- 
consistent calculations employ effective masses well be- 
low 1, the corresponding single particle energies are more 
widely spread than the experimental separation energies. 
Therefore the self-consistent approaches, in their present 
form, do not reproduce the experimental separation en- 
ergies. A precise knowledge of the separation energies 
may be very important, however, as shown in Ref.[10(. 

In TFFS the single-particle energies and wave- 
functions are the solutions of a phenomenological single- 
particle Hamiltonian and the ph-interaction is parame- 
terized universally for all nuclei. The RPA results, es- 
pecially the non-collective solutions, depend sensitively 
on the single-particle spectrum. For that reason TFFS- 
results are closer to the experimental data as, in gen- 
eral, one uses the experimental spectrum as far as it 
is available. The experimental spectrum as input into 
the RPA is of crucial importance for states that are not 
very collective, such as the odd-parity (magnetic) states. 
These states are, in general, dominated by one single ph- 
configuration and are only little shifted in energy com- 
pared to the corresponding uncorrelated ph energies. Fa- 
mous examples in this respect are the 12~ and 14~ states 
in 208 Pb [l3|]. Therefore, if one intends to do a consis- 
tent nuclear structure calculation in which collective and 
non-collective states agree with the data, one has to be- 
gin with a single-particle spectrum that is close to the 
experimental one. 

In order to investigate this problem, we will rely on an 
extended version of the Landau-Migdal theory, ETFFS. 
Landau's theory shares two important ideas with the 
modern effective field theories due to Weinberg: one has 
to identify (I) the degrees of freedom relevant for the en- 
ergy scale one is interested in and then can incorporate 
the physics of the larger energy scales in a few low-energy 
constants and (II) a small expansion parameter. How- 
ever, a systematic counting scheme to evaluate higher 
order corrections was not developed. 

In Landau's theory, the relevant degrees of freedom 
are called quasi-particles and the small parameter is the 
ratio excitation energy / Fermi energy 14]. In finite 
systems Migdal identified Landau's quasi-particles with 
the single-particle excitations in odd mass nuclei and the 
associated energy scale with the particle-hole gap, i.e. 



16 MeV (for medium mass nuclei) and 8 MeV (for the 
heavy ones), whereas the Fermi energy in nuclei is of the 
order of 40 MeV. In the TFFS only these degree of free- 
dom have been considered. On the other hand, there 
is additional degree of freedom, the phonons (l5l , [l6| . In 
the ETFFS, the phonon degree is incorporated explicitly. 
The quasi-particles are dressed by the phonons, and it is 
the dressed quasi particle energy which has to be identi- 
fied with the experimental separation energy. Details can 
be found in Ref . (l7l [Lsl [l9| (and references therein) . The 
Quasiparticle Phonon Model (QPM) by Soloviev [151 ] has 
been recently applied by Lenske et al. [2(| to the pygmy 
resonances. In contrast to the original model, here the 
authors treat the mean-field part microscopically by in- 
corporating HFB results as input for the QPM calcu- 
lations. This model allows to consider higher phonon 
states. The QRPA plus phonon coupling model by Bor- 
tignon [21] is very similar to the present one. Here the 
author start with a HF-BCS mean field were they used 
a Skyrme type interaction. In contrast to our approach 
the single particle continuum in both models is included 
in a discretized way. 

In the following we will give a derivation of the ba- 
sic equations of the TFFS as well as the extended the- 
ory within the many-body Green function method. An 
important issue is the definition of the one-body poten- 
tial from which one derives the single-particle properties, 
which are the crucial input data for the RPA. We will 
apply both approaches to the electric dipole states in 
2Q8 Pb and compare the results. We will show explicitly 
that in ETFFS the low-lying phonons give rise to a com- 
pression of the ph spectrum up to several MeV. This is 
an important result for self-consistent approaches: as the 
number of phonons and their energy are quite different in 
the various mass regions, one has to include these effects 
explicitly in self-consistent calculations. The mean field 
in self-consistent calculations provides the bare single- 
particle spectrum, which enters in the extended theory. 
Therefore, in self-consistent calculations, the lplh RPA 
is insufficient. One has to include the effects of the low- 
lying phonons explicitly, as is done in the extended the- 
ory. We demonstrate our findings by investigating the 
low-lying and high-lying electric dipole states in the lplh 
RPA as well as in the extended theory. 



II. METHOD 

The Landau Migdal theory can be derived within the 
theory of many-body Green functions [22l . [23l . I24 ]. The 
one-particle and two-particle Green functions are defined 
as 

W 2 (t 1 :t 2 ) = (-i)(A0\T{a„ 1 (t 1 )a+(t 2 )}\A0) (2.1) 



(-if (A0|T{ atfl (ti)o^(t 3 )o+(*4)a+ (* 2 )} \A0) (2.2) 
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with the time-dependent creation and annihilation oper- 
ators 



and 



a v {t) 



iHt 

e a v 



-iHt 



a+(t) = 



(2.3) 



The symbol T denotes the time-ordering operator, which 
means that the operators should be taken in time-ordered 
form with the latest time to the left and the earliest to 
the right. The nucleons are in the single-particle state 
<j) v . Here |A0) is assumed to be the exact ground state of 
an even-even nucleus. In the Landau Migdal theory one 
investigates the response function L, which is defined as: 



L(13,24) = 5(13,24) - 3(1,2)3(3,4). 



(2.4) 



= -x 



K 



FIG. 1: Graphical representation of the Bethe-Salpeter equa- 
tion for the response function in the ph channel 

The response function obeys an integral equation of 
the form 

L(13,24) = -0(1,%(3,2) 

-i J d5d6d7d8 g(l, 5)if (57, 68)i(83, 74)3(6, 2), (2.5) 

where K is an effective two-body interaction. A graphical 
representation of eq. (|2.5|) is given in Fig[T] One should 
note that the kernel K of the integral equation is the 
irreducible part of the response function with respect to 
the particle-hole propagator. We further introduce the 
vertex function f by 

L(13,24) = J d5d6 g(l, 5) r(53, 64) 3(6, 2). (2.6) 

The corresponding equation follows from eq. (|2.5p 
r(13,24) = -5(1, 4)5(3,2) 
/ / d5d6d7d8if(17, 28)3(8, 5)3(6, 7)?(53, 64). (2.7) 



In this paper we investigate the excitation energies and 
transition amplitudes of double closed-shell nuclei. With 
the time order t$,ti > ti,t% we insert into the two- 
particle Green function a complete system of eigenfunc- 
tions of the A-particle system between the two particle- 
hole pairs: 



3(13, 24) = £ 3°"(34)3 n0 (12). 



n=0 



with 

3 0,l (34) = -i(A0\T{a(3)a+('i)}\An) 



(2.8) 



(2.9) 



3™°(12) = -i(An\T{a(l)a+{2)}\A0} 



(2.10) 



The spectral representation of the response function 
eq. (|2.4p is also given by eq. (|2.8p . but the sum begins at 
n = 1 



L(13,24) = £ 3°"(34)3 n0 (12). 

n=l 



(2.11) 



With the appropriate time-ordering and a Fourier trans- 
formation one obtains the spectral representation the re- 
sponse function 



71=1 



.n0*_. ,n0 
A12 A43 



\21 A34 



^ \tt + E n - iS Q- E n + iS 



.(2.12) 



(5^0+ 



Here E n are the (exact) energies of the excited states of 
the A-particle system and 



X™ 2 ° =<An\aia 2 \AQ> 



(2.13) 



are the corresponding transition matrix elements between 
the ground state and the excited states of the A-particle 
system. In the actual calculation we solve the equation 
for the change of the density in an external field and 
determine its poles. 



A. One-particle Green function and the nuclear 
shell model 



In order to solve the equation for the response function 
and vertex function, respectively, we have to know the 
one-particle Green functions 3(1,2) and the effective two- 
body interaction K . As Landau has shown, one does not 
need the full information included in 3, but only the so- 
called pole part of the one-body Green functions near 
the Fermi surface. The rest gives rise to a renormalized 
interaction and effective one-body operators. We first 
discuss the pole part of the one-body Green function. 
For that we consider the Dyson equation in coordinate 
space [22j 



J d& £(&,6,e) 9(6,6,6 



) = 6(Zi ~ 6) , (2.14) 



where £ = (r, s) represents space and spin coordinates. 
In an arbitrary single-particle basis <p v {£) the Dyson 
equation has the form 



2m 



E(6,6,e) 



(2.15) 



We now chose a special basis <y9„(£) that diagonalizes the 
expression in the brackets 



= E Vl {e)8 v 



(2.16) 
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It is obvious that such a basis must depend on the energy 
e. The one-body Green function becomes diagonal in this 
basis and can be written as 

(2.17) 



™ e-E Ul (e)- 

Due to Landau's renormalization procedure one needs 
the one-particle Green functions only near the dominant 
poles e Vl i.e. near the poles, which are the solutions of 
the equation 



E v (e v ) 



(2.18) 



and having maximal single-particle strengths (see be- 
low). The single-particle energies e v so defined are the 
quasi-particle energies (in the sense of Landau) for finite 
systems. In the vicinity of e v we obtain for the one- 
(quasi)particle Green function 



(2.19) 



The residue of g at the pole is called the single-particle 
strength 



1 



""(l-^)e= 6 



(2.20) 



Within the framework of this method it is easy to de- 
rive an equation for the single-particle energies and the 
single-particle wave functions. To this aim let us omit, 
for the sake of simplicity, the spin variables and introduce 
the self-energy in a mixed representation: 

E(Jl, k, e) = J ds E(i2 - |s, R+ ±s, e) e isk , (2.21) 

where R = + r<z), s = r"2 — r\. The local energy 
approximation (see [391 ]) consists in the replacement of 
the function E(i?, k, e) by the following ansatz: 

E(R,k 7 e) ~ t(R, k 2 F ,e F ) 

+ (k 2 - k 2 F ) (<9E(i£, k 2 ,e F )/dk 2 ) k2=kl 

+ (e-e F )(dt(R,k%,e)/de)^ eF , (2.22) 

where E(i2, k 2 ,e) is the average with respect to the an- 
gular variables of vector k, 



E(R, k 2 ,e) 



dsE(R-ls,R+ls,e)^^- (2.23) 



and the Fermi energy e F is related to the Fermi momen- 
tum k F (r) by the formula 



2m 



k 2 F (r) + 'E(r,k 2 F (r),e F ) = e F . 



(2.24) 



On substituting eq. (|2.22[) into the Dyson equation (|2.14p 
we come, after a series of transformations, to the follow- 
ing equation for the eigenfunctions: 



2m* (r) 



V + U c (r) 



ipx(r) = e x Px(r) , 



(2.25) 



where 



m*(r) 



and 



a(r) 



2m /9E(r,fc 2 ,e F ) 
7j2" I dp 



dt(r,k 2 F (r),e) 
de 



(2.26) 



(2.27) 



If we also take the spin degree of freedom into account 
then, proceeding in the same way, one obtains the follow- 
ing single-particle Hamiltonian: 

H =- V 2^) V + U ^ 

+| [a(r) (crl) + {at) a(r)] . (2.28) 

Here I is the angular momentum operator. The form 
of this Hamiltonian follows from the expansion of the 
self-energy E about the Fermi energy e F and the Fermi 
momentum p F , and symmetry arguments. This Hamil- 
tonian agrees with phenomenological shell model Hamil- 
tonians. It is useful to point out [22\ that the spin-orbit 
term follows here from the expansion of the non-locality 
of E. A certain fraction of the non-locality may be con- 
nected with the uj and a exchange of the bare nucleon- 
nucleon interaction but certainly not the whole effect, as 
is assumed in some of the relativistic effective mean field 
approaches. 



B. Landau's renormalization procedure 

We now return to the one-particle Green function. 
With eqs. (|2.19l and l2.20"]) we can separate the one-particle 
Green function into a singular part and a remainder: 



with 



to 



+ ir) sign{e Vl — //) 



(2.29) 



^0 - E Vl 



A 

-i 



,2 



(2.30) 



z Vl = \(MS\a Vx \A+\v x )\ 
for 

-\v x )\ 

for e Ul <[i, (2.31) 

where \A ± 1 v\) denotes an excited state of the A ± 1- 
particle system. With this ansatz one writes the prod- 
uct of two one-particle Green functions — a quantity that 
appears in all the integral equations we have derived 
before — as a singular part S and the remainder i?, us- 
ing the ansatz in eq. (|2.29[) 



(2.32) 
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with 



~V3 ~ c I/4 



^1/2 ^1 ^ 



(2.33) 



Here n„ are the quasi-particle occupation numbers (1 or 
0) and f2 is the energy transfer between particle and hole 
states. Using Landau's renormalization procedure, one 
obtains from eq. (|2.8p an equation for the renormalized 
vertex function r: 



^I/l 1/3,1/21/4 O*') ^) 



V\ 1/3 ,1/21/4 



(w, n) + 



£l/ 5 Ci/g 



7"l/K I 



The renormalized vertex is defined as 



' 1/11/3,1/21/4 



(w,f2) = v / i^z 



i/ 2 '1/1^3,1/21/4 



(2.34) 



(2.35) 



Here only the singular part of the product of the two 
Green functions appears explicitly, whereas the remain- 
der gives rise to a renormalized effective two-body in- 
teraction F ph and a renormalized inhomogeneous term 
t". In similar fashion one can renormalize the response 
function L and obtains 

L(13, 24) = r w (13, 57)Z(57, 68)1^(68, 24) 
+1-^(13,57)^(57,24), (2.36) 

where L is the quasi-particle response function 

L(13,24) =S(1,2)t(1324). (2.37) 

C. The renormalized equations of the TFFS 

As we have seen in section II, the response function 
includes the transition amplitudes between the ground 
state and the excited states of an A-particle system. Us- 
ing the arguments in section II and the projection proce- 
dure described in [23j , one obtains from the renormalized 
response function in eq. (|2.36[) and the renormalized ver- 
tex function eq. (|2.34p an equation for the quasi-particle 
quasi-hole transition matrix elements 

( e "l ~~ e i/2 — ^) Xi/!i/ 2 = ( n i/l — ^2) ^ ] ^11/41/21/3 Xj/ 3 i/ 4 ; 

1/31/4 

(2.38) 

which are connected with the full ph-transition matrix 
elements X^° 2 defined in ea. (|2.13|) by the relation 



mO _ a; r, 

Al/ll/2 / j V 3 ViVil>2 



(2.39) 



The transition matrix element of a one-body operator Q 
is given by 



(Am\ Q I AO) 



- O v 

/ ^ ^1/12/2 A I/1I/2 5 
l/ll/ 2 



(2.40) 



where Q is the renormalized one-body operator 



off _ u 

V\V<2 / j 'l/i!/ 3 I/ 2 l/4 

I/3 1/4 



(2.41) 



There exist many derivations of the RPA equations 
(I2.38[) . In all cases the equations have the identical form 
but the approximations that led to them are quite dif- 
ferent. In the present derivation within the many-body 
Green functions no approximations have been made so 
far. We arrive, however, at a renormalized ph interac- 
tion which, in principle, can be calculated starting from 
the bare nucleon-nucleon interaction. This, however, is 
not done in practice. An important result of our deriva- 
tion concerns the single-particle energies, which are ma- 
jor input data. The quasi-particle energies in the single- 
particle Green functions (eq. (|2.29p ) are the experimental 
single-particle energies of the Ail particle system. In the 
case of self-consistent calculations this means that one 
has to choose an interaction that reproduces in mean field 
approximation the experimental single-particle spectrum 
of the neighboring nuclei. This point will be discussed 
again in connection with the extended theory. 

From eq. (|2.34p one obtains the change of the quasi- 
particle density in an external field (details are given in 
Ref.dlil): " 



Fy, V „ V„ V ,PV4 I/3 (O) ■ 



I/3I/4 



(2.42) 



The actual calculations presented here have been per- 
formed in r-space because this allows treatment of the 
continuum in the most efficient way. The corresponding 
equation takes the form [lj| [25| 

p(r, fi) = - J d 3 r'A{r, r', fl)Q cS {r', ft) (2.43) 

d 3 r'd 3 r"A(r, r', n)F ph (r\ r")p{r", fi) . 

The poles of the equation are the excitation energies of 
the A-particle system and p(r, fi) at a given pole is the 
corresponding transition density. 



III. EXTENDED THEORY OF FINITE FERMI 
SYSTEMS (ETFFS) 

As mentioned earlier, the original TFFS allows one 
to calculate only the centroid energies and total transi- 
tion strength of giant resonances because the approach 
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is restricted to lplh configurations. In order to describe 
nuclear structure properties in more detail one has to in- 
clude 2p2h, or even more complex configurations. A the- 
oretical approach that takes into account the complete 
2p2h configuration space and a realistic 2p2h interaction 
is numerically intractable if one uses a realistically large 
configuration space. For this reason the main approxi- 
mation in the ETFFS concerns the selection of the most 
important 2p2h configurations. One knows from exper- 
iments, e.g. from the neighboring odd mass nuclei of 
208 Pb, that the coupling of low-lying collective phonons 
to the single-particle states gives rise to a fragmentation 
of the single-particle strength, which is seen in even-even 
nuclei as a spreading width in the giant resonances. 

In the past 15 years an extension of the TFFS that 
includes, in a consistent microscopic way, the most col- 
lective low-lying phonons has been developed by some of 
us [H, QjJ [2|| (and references therein) . Many of these 
ideas have been developed in the early work by Werner 
[13 and Werner and Emrich [H (see also Refs-^HH) 
and Ref. [3(| ■ In this way one considers a special class of 
2p2h configurations because the phonons are calculated 
within the conventional (lplh) TFFS. The phonons give 
rise to a modification of the particle and hole propaga- 
tors, the ph interaction and the ground state correlations. 



(b) 



Id, 



FIG. 2: Graphs corresponding to the propagator of (a) the 
RPA and (b-d) the extended theory. The wavy lines and the 
thin lines denote the phonons and the single-particle propa- 
gators, respectively. 



quasi-particles, which now has the form 

2 



= 1 



(3.46) 



Here gf, xt (e) is the one-quasi-particle Green function in 
the extended theory and the vertex 7^f ;l , which couples 
the quasi-particle state v to the core excited configuration 
jjL®i, is given by 



7?' 



EF ph v l 

a, /3 



(3.47) 



where x 1 is the RPA wave function of the phonon con- 
sidered. The corresponding energy-dependent correction 
for the ph-interaction has the form 



A. Some basic relations 

In Landau's theory the self energy E is irreducible in 
the one-particle and one-hole channels, respectively, and 
the kernel K in the integral equation for the response 
function (see FigQJ is irreducible in the particle-hole 
channel. One now introduces a hierarchy of energy de- 
pendencies. As one can see from eq. (|2.33l) . the particle- 
hole propagator introduces a strong dependence on the 
energy transfer fi. Compared to this singular behavior, 
one neglects in Landau's theory the energy dependence of 
the ph-interaction. One neglects consistently the energy 
dependence of the self-energy E in the Dyson equation 
(|2.14[) and considers the quasi-particle and quasi-hole 
poles only. This approach is of leading order in the energy 
transfer fi. In the extended theory, in which phonons are 
introduced, one considers the next-to-leading order in the 
energy transfer; i.e. the self energy and the ph interac- 
tion become energy dependent. We then write, explicitly, 



E„(e) = E„ + E« 4 (e) ; 



pph.ext /q\ 



(3.44) 



(3.45) 



First we discuss the so-called g 2 approximation [T3|, 
which is graphically shown in FigJ5] The upper part of 
graph (b) gives a correction to the Dyson equation for 



(7^*7^ 



e - e' + (fi; - id) 



e - e' - (fij - iS) 



(3. 



This model has been applied, for example, to Ml res- 
onances in closed shell nuclei [3l[. In general, the ap- 
proximation has problems with the so-called second-order 
poles, which give rise to a distortion of the strength func- 
tion near these poles. For this reason the model has 
been extended in several steps: (I) summation of the g 
terms [32| and (II) partial summation of diagrams [lg|, 
which is called chronological decoupling of diagrams. In 
the latter case all lplh(£>phonon contributions are consis- 
tently included and all more complex configurations are 
excluded — as long as one neglects ground state correla- 
tions. The actual formulas, however, include the ground 
state correlations completely. 

In the present analytical approach, as in [18[ , two types 
of ground state correlations are included: the conven- 
tional RPA ground state correlations, which affect the 
location and the magnitude of the residua of the ph prop- 
agators only; and the new type of ground state corre- 
lations, caused by the phonons. These correlations are 
qualitatively different from the conventional RPA corre- 
lations because they create new poles in the propagator, 
which then cause transitions between the lplhgiphonon 
in the ground state and the excited states. They give rise 
to a qualitative change of the strength distribution and a 
change in the sum rules for the moments of the strength 



7 



function. In the present numerical application we only in- 
cluded the conventional RPA correlations. Calculations 
which include the complete ground state correlations are 
in progress. 

As in the TFFS, the energy dependence of the general- 
ized propagator, A ext , is much stronger than that of the 
generalized ph interaction, as the next-to-leading order 
energy dependence is removed from the interaction and 
explicitly taken into account in the generalized propaga- 
tor. Therefore one considers explicitly only the energy 
dependence of the propagator and neglects the energy 
dependence of the ph interaction. The interaction is pa- 
rameterized as before, although the corresponding pa- 
rameters may differ from the previous ones. The final 
equation for the change of the quasi-particle density in 
an external field has a form identical to that in the TFFS 
(eq.mSi): 

p ext (r, ft) = - J d 3 r'A ext (r, r', ft)Q e// (r', ft) (3.49) 

d 3 r'd 3 r"A ext {r, r' , Q)F ph '(r', r")^ xt {r' 1 ', ft) . 

The analytic form of the generalized propagator can be 
found in Ref. [H. 

In the following we investigate, within the continuum 
random phase approximation (CRPA) and the ETFFS, 
electric dipole states in 208 Pb. The strength distribu- 
tions for isovector as well as isoscalar transitions are 
shown. The single-particle continuum has been included 
at the RPA level where it is taken into account correctly 
within our Green function technique in the coordinate 
representation. For details see In our approach we 
consider therefore the three mechanisms that create the 
width of giant resonance, namely (I) the Landau damp- 
ing ((Q)RPA configurations), (II) the escape width (the 
single-particle continuum) and (III) the spreading width 
(phonon coupling, or complex configurations). As in all 
our previous calculations within the ETFFS, we include 
the most collective low-lying phonons. These phonons 
have been microscopically calculated within the RPA. In 
the present and in our previous calculations we used the 
procedure developed in [33| to obtain the energy of the 
spurious El state to be exactly equal to zero without 
having to use the procedure of fitting force parameters. 



B. Single-particle basis 

We have seen that in the conventional TFFS Landau's 
quasi-particle are the single-particle states of the neigh- 
boring odd mass nuclei. Therefore one has to use as input 
data in the corresponding equation for the excited states 
of the even nuclei (renormalized RPA) the experimental 
single-particle energies. These single-particle energies in- 
clude obviously the effect of the coupling to more compli- 
cated configurations e.g. phonons. In the actual calcula- 
tion within the extended theory we consider about twenty 
low-lying phonons explicitly and for that reason one has 
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energies [MeV] 


2+ 


4.09 6.44 


3' 


2.61 


4+ 


4.32 5.44 6.00 


5" 


3.20 3.66 5.29 6.14 7.22 


6+ 


4.42 5.21 


7- 


4.04 4.69 5.03 5.66 6.33 7.17 


8+ 


4.61 4.99 5.22 6.23 6.43 



TABLE I: Phonons used in the ETFFS calculation 



first to determine bare single-particle energies that do 
not include the coupling to those phonons. If we cou- 
ple the corresponding phonons to the bare single-particle 
states we reproduce the original quasi-particle energies. 
In order to obtain the new refined single-particle basis 
e v one has to solve eq. (|3.50|) . which can be obtain from 



E 



ft,; 



(3.50) 



(0.82) 2g 7/2 Neutron s.p. energies] 



(0.89) 3d 3/2 ' 



4s 1/2 (0.95) 



3d 5/2 (0.93) 



1j 15/2 (0.77) 



1i 11/2 (0.91) 



2g 9/2 (0.92) 



Equation ( I3.50|) has been first derived by Ring and 
Werner [3(j. It has been applied to the calculation of 
the fragmentation of the single particle strength in the 
neighboring odd mass nuclei of 20S Pb [30j and the high- 
spin states in 208 Pb Ref. [13]. 



> 




(0.94) 3p 1/2 



(0.92) 2f 5/2 



(0.92) 3p 3/2 



(0.86) 1i 13/2 



(0.78) 2f 7/2 



(0.76) 1 h 9/2 
bare energies 



exp. energies 



FIG. 3: Neutron single-particle energies in 208 Pb. Here we 
compare the experimental values with the bare single-particle 
energies which were obtained from eq. (|3.50[ ). The values in 
parentheses are the single particle strengths. 
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FIG. 4: Same as in Fig. [3] for the proton single-particle ener- 
gies in 208 Pb 



The phonons that are explicitly considered in the 
present calculation are given in Table |TJ In Figs. [3] 
and U the bare single-particle spectrum for protons and 
neutrons are shown and compared with the experimen- 
tal spectrum. The energy shift due to the coupling to 
the phonons is nearly twice as large for the protons as 
for the neutrons. In addition, one obtains a fragmenta- 
tion of the single-particle strength. The corresponding 
values are also given in Figs. [3] and |U Here we point 
out that, in addition to the phonon coupling, the short 
range correlations in the G-matrix also reduce in a uni- 
form way the single-particle strength. This part of the 
reduction is taken care of by the renormalization proce- 
dure in section II. The bare energies may be identified 
with the mean energies that one obtains within a self- 
consistent approach. From the present calculation one 
observes that the (in general) too large spreading of the 
self-consistent single-particle spectra compared with the 
experimental ones would be reduced on the average for 
the neutrons by about 1 MeV and for the protons by 
1.5 MeV with relatively large fluctuations. If one solves 
equation (|3.50|) with the corresponding mean field so- 
lutions one obtains the dressed single-particle energies 
that should agree (from our point of view) with the ex- 
perimental spectrum. The number of phonons can be 
restricted to the most collective ones because they are 
the most important ones for the width of the giant res- 
onances. They also influence the different single-particle 
states in an individual manner. If we increase the num- 
ber of phonons one only obtains an overall shift that is 



equivalent to a change in the effective mass. We conclude 
that for a given set of phonons one has to chose an effec- 
tive interaction with the appropriate effective mass and 
spin orbit interaction. We return to this question in the 
final discussion. 



C. The residual particle-hole interaction 

As in our previous calculations we used the effective 
particle-hole Landau-Migdal interaction 



F(r,r') = C [f(r) + f'(r)T 1 -r 2 
+ (.9 + 9't~i ■ T 2 )cr 1 ■ cr 2 ] 6(r - r') 



(3.51) 



with the conventional interpolation formula, for example, 
for the parameter /, 



f{r) = fex + (fin ~ fex)po(r)/po(0) 



(3.52) 



and similarly for the other r-dependent parameter f'(r). 
Here po(r) is the density distribution of the ground state 
of the nucleus under consideration and /j„ and f ex are 
the force parameters inside and outside of the nucleus. 
The standard values of the parameters, which have been 
used for all the nuclei under consideration [l9l |. are: 



fin = -0.002, f ex = -1.35, f' ex = 2.30, /' = 0.76, 



9i 



g ex =0.05, g' ir 



tie 



0.96. 



(3.53) 



The parameters are given in units of Co = 300 MeVfm 3 . 
For the nuclear density po (r) in the interpolation formula 
we chose the theoretical ground state density distribution 
of the corresponding nucleus, 



(3.54) 



Here Ri(r) are the single-particle radial wave functions 
of the Woods-Saxon single-particle model. This is more 
consistent than the previously used phenomenological 
Fermi distribution. For that reason we readjusted in the 
present calculation f ex to reproduce the breathing mode 
and obtained the value of fex = -1.35. For other details 
of the calculations, see Ref . fig. |2"1 . 



IV. LOW-LYING ELECTRIC DIPOLE 
STRENGTH 

A. The isovector case 

The giant electric dipole resonance (GDR) is one of the 
most collective states in nuclei. It exhausts a major part 
of the energy-weighted sum rule; the excitation energy 
is a smooth function of the mass number; and the form 
of the GDR changes only little from nucleus to nucleus. 
Phenomenological collective models describe well the A- 
dependence of the energy and the strength with only a 
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FIG. 5: Electromagnetic El strength distribution calculated 
within the conventional lplh RPA. In (a) we show the un- 
corrected ph-matrix elements. In (b) we used our universal 
isoscalar interaction and the spurious state has been removed. 
For the results shown in (c) a reduced isovector force has been 
used, while (d) shows the results with the full interaction. 



few parameters. In microscopic models these collective 
states are coherent superpositions of many particle-hole 
states. The theoretical results of such models for the 
mean energy and the total strength agree, in general, very 
well with the data. The strength distribution, however, 
is reproduced by only a few very involved models, such 
as the ETFFS considered here. 

Here we first investigate the question whether the 
pygmy dipole resonances (PDR) are like the GDR collec- 
tive states in the microscopic definition. That is, are they 
coherent superpositions of many particle-hole states? In 
order to answer this question we performed a series of cal- 
culations within the continuum RPA in which we varied 
the particle-hole interaction. 

In FigOi the dipole distribution for the electromag- 
netic £T-operator is plotted, where the particle-hole in- 
teraction is zero; i.e. the figure shows pure proton and 
neutron ph-matrix elements. The largest matrix ele- 
ments are the ones between the ph states with the largest 
angular momenta that differ by one unit — the stretched 



configurations. The state at E «6.5 MeV is the neu- 
tron (Iji5/2)(l«i3/2) _1 configuration and the states at 
E w 7.15 MeV and E w 7.70 MeV are the proton con- 
figurations (lii3/2)(l/in/2) _1 and (l/i 9/2 )(lg 7 /2 ) _1 . The 
three largest transitions exhaust about 40% of the to- 
tal strength. Here, we have to bear in mind that the 
transitions shown also include spurious components that 
give rise to the spurious isoscalar state at zero energy. 
This state corresponds to the translation of the whole 
nucleus. For that reason one has to remove the spurious 
components in each of the particle-hole configurations. 
The method developed in Ref. 33] is used that brings 
the spurious state exactly to zero energy and removes 
completely the spurious strength in all excited states. 

In FigJ3J>d the electromagnetic dipole strength, from 
which the spurious strength has been removed, is shown. 
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FIG. 6: Integrated electromagnetic El strength calculated 
within the ETFFS. These results should be compared with 
the distribution shown in Fig[5jl. Note that the strength in 
the present figure is integrated over the smearing width, given 
m e fm . 

To demonstrate the effect of the repulsive isovector in- 
teraction we have chosen three different values for the 
isovector force parameter f' ex . The isoscalar ph interac- 
tion is the same for all three cases. It reproduces the 
breathing mode, which is later shown in FigJ2Dl If we 
compare the completely uncorrelated strength distribu- 
tion in FigEK with the results in FiglSb, where the spuri- 
ous components are removed but the isovector force pa- 
rameter f' ex = 0, one observes an appreciable shift of the 
B(E1) strength to higher energies. We notice thereby 
that if one removes the spurious strength one has to in- 
troduce correlations that give rise to this redistribution 
of the strength. With setting the isovector interaction 
to zero, most of the B(E1) strength (87%) lies below 10 
MeV. If we increase the isovector interaction parameter 
t° f'ex = 1-0 (FigHt), which is about half of the force 
parameter that reproduces the experimental data of the 
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GDR, the dipole strength below 10 MeV is reduced to 
11%. 
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In FigtZ] the electromagnetic El strength distribution 
calculated within the ETFFS is shown up to 20 MeV with 
a smearing parameter A = 250 keV and three different 
f' ex parameters. It is interesting to see how the strength 
is shifted to higher energies with increasing interaction 
strength. For our conventional value f' ex = 2.30 the the- 
oretical distribution in the GDR region agrees quantita- 
tively with the experimental data. In connection with 
the previous discussion, the energy range between 6-8 
MeV is of special interest. With no isovector interac- 
tion, but the spurious components removed, one obtains 
some concentration of strength around 6.2 MeV and 7.5 
MeV. The latter is considered as the PDR region. With 
increasing isovector force this strength is reduced. For 
the full force the reduction is a factor of 5 compared to 
f' ex — 0. This behavior is just the opposite of a collective 
structure, where with increasing force the collectivity is 
enhanced. 



FIG. 7: Electromagnetic El strength from 5 to 20 MeV cal- 
culated in our extended theory. The thick line indicates the 
final result whereas the other lines are the results of calcula- 
tions where the isovector force parameter has been varied as 
indicated in the figure. 

Finally, in FigfSJi, where we use our conventional force 
parameter (which reproduces nearly quantitatively the 
experimental data of the GDR), only 5.2% of the strength 
remains below 10 MeV. We observe that the major part 
of the El strength is shifted into the GDR region, where 
it creates a very collective resonance — in agreement with 
the data. 

Within our model, in the electromagnetic El strength 
up to 9 MeV, which is the region of the pygmy reso- 
nance, we do not find one single state in which most of 
the low-lying strength is concentrated, but we find al- 
ready in the lplh RPA several states whose energies are 
essentially unchanged compared with the f' ex = case. 
The El strength, however, is reduced by a factor of 10 
due to the strongly repulsive isovector force, which shifts 
the strength into the GDR region and produces a col- 
lective resonance. Such a behavior is expected from the 
schematic model of Brown and Bolsterli [36[ . 

In FigJS] we present the results of the ETFFS calcula- 
tion where, in addition to the single-particle continuum, 
the effect of the phonons is also included. The phonons 
that we consider are given in Table HI The result in Fig (HI 
should be compared with the results in FigJSJi. One 
observes a further fragmentation of the strength and a 
small shift of the two strongest states to higher energies. 
Whereas the calculated strength distribution between 7 
and 8 MeV is in good agreement with the data, we obtain 
too little strength below 6 MeV and too much between 
6-7 MeV as compared with the present data. A shift of 
the two neutron states I715/2 and (1*13/2) 1 improves the 
agreement between theory and experiment. 
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FIG. 8: Summed transition densities for the electric dipole 
response shown in Fig[6]from 5 to 6 MeV 

In Figs l8lfT0l the transition densities between 5-6 MeV, 
6-7 MeV and 7-8 MeV of the electric dipole response are 
plotted for transitions shown in Fig[6] In all three cases 
the protons and neutrons are in phase, which indicates 
strong isoscalar components. This at-first-glance surpris- 
ing behavior has also been found by other authors. It can 
be understood immediately in our microscopic model. As 
discussed in the next section, the isoscalar force is at- 
tractive and gives rise to a collective structure around 22 
MeV. Due to this attractive force, some of the high-lying 
strength couples also to the low-lying states, whereas 
the strong repulsive isovector force simultaneously shifts 
the isovector strength from these configurations into the 
GDR region. We demonstrate this behavior in Figs fT2lfT^l 
where the transition densities are plotted for the same en- 
ergy ranges but with the force parameters f set to zero. 
As the isovector force is zero, the isovector strength re- 
mains in the low energy region, which can be clearly seen 
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FIG. 9: Summed transition densities for the electric dipole 
response shown in Fig[6]from 6 to 7 MeV 



FIG. 11: El isoscalar strength distribution in Pb from 5-10 
MeV for the transition operator (r 3 — 5/3(r 2 )r)Yi iM . Here we 
compare the results derived from the CRPA with the results 
obtained with the extended theory. 
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FIG. 10: Summed transition densities for the electric dipole 
response shown in Fig[S]from 7 to 8 MeV 

in the transition densities. 



B. The isoscalar case 

As mentioned earlier, for the electric isoscalar dipole 
states we have a special situation because the lowest 
isoscalar resonance is the spurious state corresponding 
to the translation of the whole nucleus. In self-consistent 
calculations this state appears (at least in principle) at 
zero energy and carries all the spurious strength. 

In the Landau-Migdal approaches the spurious state 
has to be removed explicitly. Here we use the procedure 
developed in Ref.[33| to obtain the energy of the spu- 



FIG. 12: Summed transition densities from 5 to 6 MeV with 
the isovector force parameter f' ex — 



rious El state exactly at zero energy, with no spurious 
components in the excited state. As the isoscalar electric 
dipole strength distribution depends on the isoscalar pa- 
rameter /, we test the force by calculating the isoscalar 
monopole resonance — the breathing mode. With our 
universal value of f ex = —1.35 we reproduce the mean 
energy as well as the width of the resonance (see the 
next section). In FigQT] the results of our calculation 
for the isoscalar dipole operator (r 3 — 5/3(r 2 )r)Yi, p up 
to 10 MeV are presented. There we compare the CRPA 
with the ETFFS results. In the CRPA we obtain one- 
strong state near 7.4 MeV and several somewhat weaker 
states at lower energies. Due to the phonon coupling the 
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FIG. 13: Summed transition densities from 6 to 7 MeV with 
the isovector force parameter f ex — 
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FIG. 14: Summed transition densities from 7 to 8 MeV with 
the isovector force parameter f' ex — 



strength of the strongest state is fragmented and shifted 
to slightly higher energies. Those isoscalar states have 
been investigated experimentally with the (a, a', 70) re- 
action. So far the experimental data from 208 Pb are not 
given in the form of cross sections but only in counts 
[6|. The data show a very strong signal at about 5.6 
MeV and a somewhat weaker signal around 6.7 MeV. As 
the isoscalar strength distribution in Fig[ll]is calculated 
with the operator (r 3 — 5/3(r 2 )r)Yi iM , a comparison with 
these data is not directly possible because the (a, a') re- 
action is sensitive to the tails of the isoscalar distribu- 
tions and the 70 is proportional the isovector admixture. 
Unfortunately, analyzing programs in which microscopic 
transition densities can be used do not yet exist [jj . The 



transition densities shown previously are the same for the 
isoscalar and isovector case. Only the transition strength 
is different because the corresponding operators weight 
the transition densities differently. 



HIGH-LYING ELECTRIC DIPOLE 
STRENGTH 
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FIG. 15: El isovector strength from 5 to 20 MeV. We com- 
pare here the results of the continuum RPA with the extended 
theory the width and the mean energy of the latter agree 
quantitatively with the experimental results. 
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FIG. 16: Transition density of the isovector giant dipole res- 
onance 

The GDR was calculated with our universal parame- 
ters given in eq. (|3.53[) and the result is shown in FigfTSl 
The mean energy and the width of the resonance agree 
with the data within the error bars. The escape and 
spreading widths are included in our microscopic model. 
In addition, we consider a smearing parameter A = 250 
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FIG. 19: Isoscalar El strength calculated with a smearing 
parameter A= 3 MeV. Note: The larger smearing parameter 
is justified only for the high-lying (>15 MeV) part of the 
resonance. 



FIG. 17: Isoscalar El strength from 5 to 30 MeV calculated 
with a smearing parameter A = 250 keV 
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FIG. 18: Transition density of the isoscalar giant dipole res- 
onance 



keV , which corrects for more complex configurations that 
have not be considered in our approach. Here we point 
out again that the low- lying and high- lying El strength is 
calculated within the same model with exactly the same 
single-particle energies and force parameters. This differs 
from some self-consistent approaches in which the single- 
particle spectrum is modified in order to fit the low-lying 
spectrum e.g. (37j . 

In Fig ll6l the corresponding transition density is plot- 
ted. It is surface-peaked and protons and neutrons are 
out of phase. As we have 50% more neutrons than pro- 
tons, the neutron density has not only a stronger and 
longer tail than the proton density, but it is also peaked 
further out. For that reason we also have some admix- 
ture of isoscalar El strength — even in the isovector giant 
dipole region. The form of the transition density agrees 
with surface-peaked phenomenological models. 

In Fig[T7]the theoretical isoscalar El response from 5- 



FIG. 20: Isoscalar E0 strength from 5 to 20 MeV calculated 
in the continuum RPA and the extended theory. The mean 
energy and width of the latter agrees with the data. 



30 MeV is plotted. One realizes that, in addition to the 
low-lying strength, we obtain a strong collective struc- 
ture above 20 MeV, which represents the isoscalar electric 
dipole resonance. The mean energy of the resonance is 
E t h = 22.1 MeV and the width T th = 3.8 MeV if we chose 
a smearing parameter A — 250 keV. The deep-lying holes 
contribute a major part to the high-lying El-strength. 
These holes are several MeV broad. These widths are 
not included in the calculation shown in Fig |17l where 
we used a smearing parameter of A = 250 keV. In order 
to correct for these widths we performed a calculation 
with A = 3 MeV, the corresponding result is shown in 
FigHH] The theoretical results can be compared with a 
recent (a, a') experiment [Hj], which obtained E exp = 
22.5-23.0 MeV and a width T exv sa 10 MeV. The corre- 
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FIG. 21: Transition density of the breathing mode 



the experimental single-particle energies, then with in- 
creasing isovector force the isovector dipole strength is 
shifted to higher energies and one is finally left with a 
small fraction of the sum rule in a few states that are 
only slightly shifted with respect to the unperturbed ph 
energies. There is no indication in the present model of 
a low-lying collective state that is built up by the isovec- 
tor interaction. The theoretical electromagnetic strength 
below 8 MeV is 1.94 [e 2 /m 2 ] compared with the experi- 
mental strength of 1.32 [e 2 fm 2 ]. The distribution of the 
theoretical strength agrees in the range between 7-8 MeV 
with the current data. The calculated strength between 
6-7 MeV, however, is too large and the strength below 6 
MeV too small compared with the experimental values. 
On the other hand, the width and the mean energy of 
the GDR is reproduced quantitatively in our model. 



sponding transition density is given in Fig[l8j It has the 
form of a compression mode and for that reason we have 
shown — finally — in Figsl2l71andl2"l"lthe isoscalar monopole 
strength and the corresponding transition density. These 
calculations are performed for purpose of consistency be- 
cause the isoscalar monopole and isoscalar dipole depend 
both on the isoscalar parameters / and have, in addition, 
a very similar structure, as can be seen in the transition 
densities. In the monopole case the radial integral over 
the transition density has to be zero; in the electric dipole 
case the transition density multiplied by the radial coor- 
dinate r must vanish. 



VI. DISCUSSION 

The center of our investigations has been the origin 
and the structure of the pygmy resonances. We have 
chosen in 208 Pb a nucleus with an excess of 44 neutrons, 
which means that we have 50% more neutrons than pro- 
tons. The structure of this doubly-magic nucleus, as well 
as of the neighboring odd mass nuclei, is well known 
and all microscopic models and microscopic theories — 
self-consistent ones and Landau-Migdal — work best for 
this case. 

Our calculation shows, first of all, that we obtain 
isovector as well as isoscalar electric dipole states be- 
low 10 MeV. This is obvious, as the ph-interaction has 
isoscalar as well as isovector components. Here we point 
out that all the spurious components are removed; i.e. 
the expectation value of the translation operator rYi M in 
all the excited states is negligibly small. The isoscalar 
strength is due to the 3fku and higher ph components. 
In the present Landau-Migdal approach one is able to 
remove the spurious strength without performing a full 
RPA calculation, as must be done in self-consistent ap- 
proaches. Here, the force parameters can be varied inde- 
pendently This allows one to investigate the influence of 
the isovector interaction on the low-lying isovector dipole 
strength in detail. 

We have shown in Figs. 5 and 7 that if one starts with 



It is obvious that if one changes the single-particle 
spectrum one is always able to improve the agreement 
with the data. In the present calculation we only con- 
sider the lowest order in the expansion of the Migdal in- 
teraction, as is done in nearly all calculations. The next- 
to-leading order in Migdal's interaction depends on the 
momenta through pp', with the isoscalar force parame- 
ter fx and the isovector force parameter f[ . The sign and 
magnitude of the two parameters can be estimated from 
the effective mass of the nucleus and the orbital part of 
the effective magnetic operator [23| . From this consider- 
ation it follows that both parameters are negative, with 
the /j parameter slightly smaller in magnitude. These 
parameters would give rise to an attractive velocity de- 
pendent isoscalar and isovector ph interaction. A small 
effective mass, as is used in relativistic RPA calculations, 
would increase the magnitude of both parameters. A 
strong isovector, velocity dependent force could be the 
explanation of the collective states with large vorticities 
that are found in relativistic RPA calculations [ll|. In 
this approach two phonon states will give rise to a frag- 
mentation that should explain the El data between 5 
and 9 MeV. Unfortunately, the lowest two phonon states 
are at 6.7 and 6.9 MeV, whereas the lowest experimental 
El-states are more then one MeV lower. We therefore 
consider the present result, where the El strength be- 
tween 5.5-9 MeV lies essentially at the unperturbed ph- 
energies and the major part of the isovector El strength 
is shifted into the GDR, to be the more physical one. 
Moreover, a collective phenomenon has not been found 
in self-consistent calculations in which one starts with ef- 
fective forces of the Skyrme type. The next-to-leading 
order in Migdal's interaction (fi,f[) will lead to some 
redistribution of the low-lying strength which may im- 
prove the agreement with the data below 7 MeV. Such 
calculations are in progress. 
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VII. CONSEQUENCES FOR 
SELF-CONSISTENT CALCULATIONS 

The present calculation is not self-consistent in the 
sense that one starts with an effective Lagrangian (or 
Hamiltonian) , the numerous parameters of which are, in 
general, adjusted to gross properties of nuclei and can be 
used for all nuclei. These parameters cannot be deter- 
mined in a unique way and therefore there exist numer- 
ous sets of parameters that reproduce the gross proper- 
ties equally well, but give quite different results for ex- 
cited states. If we consider, as a specific example, ef- 
fective forces of the Skyrme type, we find parameteri- 
zations that, for example, give rise to effective masses 
from m*/m = .6 to m*/m = 1. For nuclear structure 
calculations the effective mass is of special importance 
because the spacing of the single-particle spectrum is in- 
versely proportional to it; a smaller m*/m expands the 
single-particle spectrum and a larger one compresses the 
spectrum. The experimental single-particle spectra in 
the neighboring odd mass nuclei of 208 Pb correspond to 
an effective mass of roughly one, whereas for medium 
mass nuclei it is smaller then one. If one describes the 
spectrum in the odd mass nuclei correctly one simultane- 
ously reproduces the non-collective states in 208 Pb. The 
collective states, on the other hand, depend not only on 
the ph-spectrum but also on the residual ph interaction. 
It seems always possible to find, under the numerous sets 
of Skyrme parameters, one set that reproduces the col- 
lective states in which one is interested, whereas the non- 
collective ones deviate from the data by several MeV. 

In our calculations we investigated simultaneously the 
El spectrum between 5-8 MeV and the high-lying spec- 
trum. Our calculations show that the El spectrum below 
8 MeV is very similar to the unperturbed experimental 
ph spectrum, whereas the high-lying collective states de- 
pends sensitively on the (universal) Migdal parameters. 
This is not an exception; rather, it is the general expe- 
rience. In the low energy spectrum of 208 Pb one finds 
some collective states with natural parity that are col- 
lective (e.g. 3~,5~,2 + etc.) but the majority of the 
states — especially the low-lying unnatural ones — are en- 
ergetically at the experimental ph energies. Therefore, if 
one is looking for a unified description of the structure of 
even-even nuclei, one has also to consider the spectra of 
the neighboring odd mass nuclei. 



VIII. SUMMARY 

The extended theory of finite Fermi systems has been 
applied to the low-lying and high-lying El spectrum of 
208 Pb as an example of a neutron-rich nucleus. In the 
present approach the low-lying electromagnetic El spec- 
trum is non-collective in the sense that the major part of 
the strength is not concentrated in one single state. The 
strength is rather distributed over several states which 
are close to unperturbed ph energies and we find in the 
low-lying El states strong isoscalar admixtures. The the- 
oretical results for the GDR as well as the breathing mode 
agree quantitatively with the data. We also find in our 
calculation a well-localized isoscalar El resonance at 22.2 
MeV, with a width of 9.3 MeV. This resonance has been 
detected in (a, a') scattering. 

From the present investigation one may draw impor- 
tant conclusions concerning self-consistent calculations. 
We have seen that the inclusion of the phonons gives 
rise to a shift of the single-particle spectrum of about 
1 MeV for neutrons and nearly 2 MeV for protons. In 
the present extended theory one has to determine bare 
single-particle energies that give the experimental ones 
if one includes the phonons. Self-consistent approaches 
determine the bare spectrum that we need to use as in- 
put in our extended theory. One therefore has to chose 
a paramctrization of the effective intercation that repro- 
duces — if one includes the phonons — the experimental 
spectrum. From this point of view, self-consistent calcu- 
lations for excited states at the level of the lplh RPA are 
not appropriate. 
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